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Abstract 

We compute the order ag QCD corrections to the density matrix for the 
decay of a neutral Higgs boson tp with (in) definite CP parity into a quark 
antiquark pair, respectively the QED corrections for the decay into a pair of 
charged leptons. We classify and calculate single spin asymmetries and spin- 
spin correlations which are generated by the scalar and pseudoscalar Yukawa 
couplings. These spin effects can be traced in tp ^ and, for heavy Higgs 

bosons, in (p —> tt. We also calculate resulting correlations among the final 
states and estimate, for the respective decay modes, the number of events 
needed to measure the Yukawa couplings with these correlations at the 3cr 
level. 
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I. INTRODUCTION 


Many extensions of the Standard Model (SM) involve more than one scalar held mnltiplet 
and thns predict the existence of more than jnst one Higgs boson. In particnlar, qnite a 
number of these models - including the 2 Higgs doublet extensions of the SM |I|-^ - allow 
for the violation of CP symmetry by the scalar self interactions. This type of CP violation is 
of great interest for scenarios that attempt to explain the baryon asymmetry of the universe 

i- 


A manifestation of CP violation in the scalar sector would be the existence of neutral 
Higgs bosons of undehned CP parity, i.e., Higgs bosons having both scalar and pseu¬ 
doscalar couplings to quarks and leptons. If neutral Higgs boson(s) should be discovered in 
the future, one would eventually like to know its (their) P and CP quantum numbers. In 
this context several proposals and theoretical studies have been made in the literature. For 
instance if (p has both scalar and pseudoscalar Yukawa couplings then a CP-violating spin- 
spin correlation is induced already at tree level in the decays of (p to fermion-antifermion 
pairs. This spin-spin correlation could be traced in (p —> and (p —>■ ft ^j. (For re¬ 

lated proposals, see |l7|-|I^.) The modes tp W^W~, ZZ could also be employed to infer 
the parity and CP properties |l^ of ip iJn||T|-|Tl. Other reactions that may be used for 
this purpose are associated ttip production ip production by high energetic photon 

photon collisions ||TI| , ^, |^ , and by high energetic annihilation p3|,^ . 


In this article we calculate the order as QCD corrections to the density matrix for 
the decay of a neutral Higgs boson ip with arbitrary scalar and pseudoscalar Yukawa cou¬ 
plings into a quark antiquark pair. (QCD corrections to the decay width of a scalar and/or 
pseudoscalar ip ^ qq were computed in [p5| -|30|.) We classify and calculate single spin asym¬ 
metries and spin-spin correlations which are generated by the scalar and/or pseudoscalar 
Yukawa couplings. These spin effects can be traced for heavy Higgs bosons in ip ^ ti 
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because top quarks decay before they hadronize and because they auto-analyze their spins 
through their parity-violating weak decays. We show that there are two correlations among 
the hnal states from ti decay whose combined use allows to investigate the CP property of 
a neutral ip. 

We then apply these results to the decay ip —*• t~t~^ and its QED corrections. We 
analyze correlations for r decay modes that have the best r spin analyzer quality. Finally 
we estimate, for the respective p decay modes, the number of events needed to measure the 
top and T Yukawa couplings with these correlations at the 3 ct level. 

II. THE DECAY ffX 

Let us briefly recapitulate the salient features of the simplest models that predict neutral 
Higgs bosons with undehned CP parity; these are 2 Higgs doublet extensions of the SM with 
natural flavor conservation at the tree level (see, e.g., [Q). Explicit CP violation in the Higgs 
potential leads to three physical neutral boson states with scalar and pseudoscalar couplings 
to fermions. In the following p denotes one of these bosons. The Yukawa interactions with 
a quark or lepton held / read: 

Cy = -^f{af + i-f5af)fp , (H.l) 

V 

where v = and Gp is Fermi’s constant, mj is the fermion mass, aj and af 

are scalar and pseudoscalar coupling constants, respectively. If a/h/ 7 ^ 0 then Cy is CP- 
violating. 

For the reaction p —> ffX the spin density matrix of the // subsystem is dehned in the 
// center of mass system by: 

/3i/32(k) = '^{f{h,ai)J{k2,f3i),X\T\p{q)) 

X 

{f{h,a2)J{k2,(32),X\T\p{q)y . (H.2) 
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Here k = ki is the / momentum in the // c.m. frame and a and (3 are spin indices. The 


sum in ([11.2|) is taken over all discrete and continuous degrees of freedom of the unobserved 
part X of the final state. The matrix ( |11.2|) can be decomposed in the spin spaces of / and 
/ as follows: 

i? = H II ® 11 + Bf cr*(8) II + B~ 11 ® cr* + Cij cr*® crh (H-S) 


The hrst (second) factor in the tensor products of the 2x2 unit matrix 11 and of the Pauli 
matrices cr* refers to the /(/) spin space. Because of rotational invariance, the functions 
Bf and Cij can be further decomposed: 

Bf = b±ki , 

Cij = Ci6ij + C2kikj + c^eijiki , (II-4) 


with k = k/|k| and coefficients A, b± and q, which only depend on the masses and couplings. 
These coefficients can all be separately measured by using suitable observables built from 
the spin operators of / and /. The trace of R, Tr{i?} = 4A is a Lorentz scalar, and is 
related to the decay rate r((p —> f fX) by: 


T = 


2m^ Svr 


Tr{fl} , 


(11.5) 


where /3 = (1 — Amj/rnffC jg velocity of the fermion in the (p rest frame. One may 
construct a complete set of observables that determine the other coefficients: 


= k ■ (si + S2) , 
(Pi = k ■ (si — S2) , 
(P2 = k ■ (si X S2) , 


(P3 — Si ■ S2 , 

04 = (k-Si)(k-S 2 ) , (II. 6 ) 
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where Si = ^cr 0 11 and S 2 = |ll ® cr are the spin operators of / and /. The expectation 
values of the observables are given by: 


{O,) 


Tr{R ■ O,} 
Tr{R} 


(11.7) 


These expectation values are trivially related to the coefficients of the density matrix (|II.2|) : 


(Oo) 

(Oi) 

(O 2 } 

(O 3 } 

(O,) 


2 ( 6 ++ 6 _) 
4A 

2 ( 6 + - 6 -) 
4A 
2c3 
T4 ’ 

3ci + C 2 

4A ’ 

Cl + C2 

4A ■ 


( 11 . 8 ) 


In Table 1 we give a list of the transformation properties of the coefficients under discrete 
symmetries. The T and CPT transformation properties hold up to non-hermitian contribu¬ 
tions to the decay amplitude. This table tells us that the CP-invariance relation 6 + = 6 _ 
can only be violated if CP is not conserved (i.e., afdj 7 ^ 0) and if absorptive parts of the 
amplitude are taken into account. In fact, the functions 6 + are zero to leading order for 
model (ini). Since the Born decay matrix is hermitian, we get from CPT invariance the 
relation 6 + = 6 ^;. Further the interaction (|II.1|) is C-conserving, which implies b± = - 6 + 
and hence 6 + = 0. On the other hand the CP-invariance relation C 3 = 0 can be - according 
to Table 1 - violated and is, in fact, violated by (|II.1|) already at tree level. 

We will now discuss some general properties of the observables (|1I.6|) . Since the ob¬ 
servable Oo is C-odd, its expectation value is zero at tree level and at arbitrary order in 
the QCD or QED couplings. CP-violating absorptive parts generated by one and higher 
loop corrections induce a nonzero difference 6 + — 6 _ and therefore a nonzero expectation 
value of the observable Oi. The CP-violating triple correlation (O 2 } is generated by ( [11. Ij) 
already at tree level and receives at the next order a contribution from the CP-violating 
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dispersive part of the one-loop decay amplitude. The other two observables 03^4 measure 
CP-conserving spin-spin correlations. 

We will now discuss the decay of a Higgs boson into top quark pairs, f = t (Later we 
will also consider the case f = t.). We have calculated the density matrix R to order 
We use the following notation: 

R = Rq H-) 

r = ro + ^r4 + o(4), 

{ 0 .> = { 0 .> +^{ 0 .) + 0 (^) , ( 11 , 9 ) 

0 Tj- 1 


where 


(o^) = 


Tr{flo ■ O,} 
Tr{Ro} 


^ Tr{R,.0.} _ ^ 

^ Tr{Ro} ^ *^oro ■ 

At Born level the rate and the expectation values are found to be: 

Nc ml 


(II.IO) 


rn = 


Stt V‘ 


-m^l3{al[3^ -\- a 


{0,}„ 

{0,\ 


0 , 

-atdtl3 

+ ~al ’ 

— 3a^ 

4(a?/?2 + ~al) ’ 
1 

“4 ’ 


( 11 . 11 ) 


where Nc = 3. At tree level the ft-system is in a pure state. This can also be explicitly 


checked using (|II.11|) . Radiative corrections to the coefficients will in general lead to mixed 
states. 

The calculation of the QCD corrections to the results ( |11.11|) proceeds as follows: We 
calculate separately the virtual and real corrections to order All singularities are treated 
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by dimensional regularization. Soft (IR) singularities appear both in the virtual and real 
corrections; they cancel explicitly in the sum of the two in accordance with the Kinoshita- 
Lee-Nauenberg theorem. We note that the necessary integrations over the hard gluon mo¬ 
mentum can be done analytically for all coefficients of the density matrix. In order to write 
the results for hi and Tr{i?i ■ Oi} entering (|II.9|) and (|II.10|) in a compact form, we define 
a; = (1 — /3)/(l + /3). Further Cf = 4/3, and Li 2 (x) denotes the dilogarithmic integral. In 
the results below the top quark mass is defined in the on-shell scheme. For the order 
contribution to the decay rate we find: 

Ti = I (a 2/?2 + 52 ^ (- 2 ^ + ln(o;) (l -F /J^)) + in(l - o;)"! 


An 

Li2(cu2) 


Li2(a;) -1 


+ (l+/9') 

,0 /i _ o 2 \ + 6/3^ -F 1) ln(a;) + 2(5{1 + (5“^) 

P } 32^2 


(13/3n48/?3-34;52 3) - 6/? (1 - 7/?^) 


32;52 


( 11 . 12 ) 


The result (|II.12|) reduces to the standard model result for = 1, Oj = 0 and agrees with 


the one given in |^. In Fig.la-c we plot the decay rate cp ^ tt for different values of at 
and dt as a function of the Higgs mass for rrit = 175 GeV. We include the 1-loop running 
of as with five active flavours and Aqcd = 200 MeV. In all plots the dashed line represents 
the Born result and the full line is the result to order Fig la and Fig Ic depict the well 
known cases of a scalar and pseudoscalar Higgs boson, respectively. For illustrating the 
case of a Higgs boson with indefinite CP parity we choose at = (2/3)^^^ and dt = (1/3)^/^. 
In the two Higgs doublet model this corresponds to vi = V 2 and maximal mixing of the 
three neutral Higgs bosons with definite CP parity (see, e.g. |^). For a given Higgs mass a 
nonzero pseudoscalar coupling enhances the rate. Note that in general the QCD corrections 
are important. 

We now turn to the other observables defined in ( [H.(j| ). The QCD corrections to their 
expectation values are determined by (cf (III.IOI) ): 
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( 11 , 13 ) 


m 


Tr{i?i ■ Oi} = -^NcCp^TTm^Jl - (3^)atat , 


mt 


Tr{i?i ■ O 2 } = ^NcCpm^^atat |2 (1 + (3^) [Li 2 (a;") - 4Li2(a;)J + 8/31n(l - a; 
+2 ln(a;) — 2f3 ln(l + 00 ) — 2f3 (l + 2/?^^ — ln^(a;) (1 — /?)^ 

+ ln(o;) (1 - 2^ - 3/?^) | , 


Tr-rp _raiNcCFml 

■ ^ 3 } - 


I ial(3‘^ - 48(1 + (3^) ( 112 ( 0 ;^) + 2Li2(a;)) 

—96(1 — /9^)Li2(a;^) + 48/?^ (^ln(a;)(l + /9^) — 2/3^ {2 ln(l — a;) + ln(l + a;)^ 
+8(1 - /32) (^3 In(o^) (In(o^) - 4 ln(l + a;)) - + 6 (23 + (3 

+3 ln(a;) (7(1 - (3^ - 24/3^ + 48/3^ + 16)' 


++ 


48(1 


2 Li 2 (ci;) — Li 2 (a;^) — 2 lrL(a;) ln(l + a;)) 


+6/3 (69 - 4/3^ - 9/3^) + 8 (vr^ - 3 \n^{uj 
+3 ln(o;) (69 - 9/3® - 43/3^ + 15/3^ 


2\2 


-4) 


Tr{7?i ■ O 4 } = /^')(Li2(o;') + 4Li2(o;)) - 48Li2(o;2) 

+48/3 (1 + /32)( ln(o; + 1) + 2 ln(l - cu)) + 12 ln2(o;) (1 - /3f (1 + /3f 
-24 ln(o;) (2(/3^ + 2/3^ + 1) ln(l - oj) + (s + (3^) ln(w + 1)) - 30/3 (l + [3'^) 
—47r^(l — /3^)(l + 2/3^) + 3 ln(a;) (l0/3^ — /3"^ + 15 — 24/3(1 + /3^)) (a^+ a^) 
+24 (1 — /3^) (1 + /3^)(Li2(ctJ^) + 2Li2(ctj)) + (ln(ci; + 1) + 2 ln(l — w)) 
x(ln(w)(l+ /3^) - 2 / 3 ) -ln(o;)(l -^^) + 3/3(l + ln(o;))](ai - a^) | . 


) 

(11.14) 


(11.15) 


(11.16) 


The strikingly simple result obtained from (p!I.13|) for {Oi) is plotted in Fig. 2, again for 
at = (2/3)^/^ and at = (1/3)^^^. Remember that a nonzero value for this correlation signals 
CP violation induced by absorptive parts. Close to the ti threshold the effect is of the order 

















of 20%. Note that close to threshold the infrared singularities due to Coulomb gluons should 
be resummed and hence our result will be substantially modihed by higher order corrections 
which we do not consider here. For Higgs masses around 400 GeV the correlation drops to 
about 10%. To exhibit the dependence of the expectation values on the unknown model 
parameters a* and dt we assume for the moment atdt > 0 and dehne 


n = 


at 


(11.17) 


at + ttt 

which takes values 0 < rt < 1. It is easy to see that the expectation values of the observables 
([II. 6 |) depend only on rt and not separately on at and dt. In Fig. 3 this dependence is shown 
for ( Oi) and a Higgs mass = 400 GeV. Around r* = 0.3 the value of (Oi) is about 
11%. The case atdt < 0 may be analysed in complete analogy by dehning rt = dt/{dt — at). 
Glearly, both cases can be distinguished by measuring the sign of {Oi). 

Fig. 4 shows the GP-violating triple correlation ( O 2 ) dehned in ( jlf.Oj) for the same set 
of parameters as used in Fig. 2. This correlation is nonzero at tree level and the QGD 
corrections are tiny. In contrast to ((Pi), h vanishes for m^—>- 2mt. For a given Higgs mass 
the correlation may reach values of ±0.5. The dependence of {O 2 ) on is shown in Fig. 
5, again for m^p = 400 GeV. We note that for all rt the order as corrections are again very 
small. 

We now discuss the results for the GP-even observables. For a scalar (pseudoscalar) 
Higgs particle, the expectation value ((P 3 ) is 1/4 (—3/4) at tree level. QGD corrections 
to these results are smaller than 1% for all Higgs masses below 500 GeV and vanish for 
rrip, = 2mt. The expectation value ((P 3 ) is quite sensitive to the value of rt, as shown in 

Fig. 6 . 

Finally, the QGD corrections to the expectation value ((P 4 ), which is —1/4 for all values 
of at and dt, are at most 2% for < 500 GeV. The dependence of these corrections on 
rt is very weak as demonstrated in Fig. 7. This completes our discussion of the order 
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corrections to the ti density matrix (| 11 . 2 |) . 

For a Higgs boson with a mass smaller than 2mw the main decay modes are ip bb, cc 
and if —>■ Becanse of hadronization effects, a spin analysis is difficnlt for b and c 

qnark pairs ~ bnt it is feasible for ip —>■ We illnstrate this below for the case of a 

Higgs boson with = 100 GeV decaying into a pair of r leptons. The QED resnlt (again 
in the on-shell scheme) follows from the formnlae ([H.9|) - (P.17|) by replacing 


N, 


c 


C, 


Q.Q 




( 11 . 18 ) 


The partial decay rate r((p —> t~t^) to order a is 0.202 MeV, which is 2.4% smaller than 
the Born rate. This resnlt is practically the same for all the three choices of the model 
parameters a,- and a,- which were also nsed before in the case oi p ^ tt. Becanse of the 
factors m^/rn^ ~ 0 ( 10 “^) and a{m^) ~ 1/128, the correlation ( Oi) is tiny and is therefore 
not a usefnl tool for analyzing the CP natnre of a light Higgs boson. In contrast, the triple 
correlation {O 2 ) takes the large valne —0.47 for Ur = (2/3)^/^ and dr = (1/3)^/^. Its 
dependence on Vr is shown in Fig. 8 . The QED corrections to this resnlt are below 1%. As 
in the QCD case, the CP-even spin-spin correlation ((P 3 ) is also sensitive to the valne 
of Tr- This is shown in Fig. 9. Finally, the QED corrections to the valne { 04 ) = —1/4 are 
below 0.3%. 


III. CP ASYMMETRIES IN ANGULAR CORRELATIONS OF THE FINAL 

STATES 

In this section we investigate the prospects to determine the CP parity of the Higgs 
boson decaying into tt or pairs by measnring snitable angnlar correlations between 

the final state particles into which the fermions decay. As we have shown in the last section, 
the radiative corrections to normalized observables are of the order of a few percent. From 
now on we will therefore neglect these corrections and perform a leading order analysis of 
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the CP asymmetries. In this approximation the rest frame of the Higgs boson is identical 
to the // c.m. system. In other words, the relation kj = —ky holds and will be used in 
constructing the CP asymmetries below. It is straightforward to generalize the correlations 
£i,£ 2 )£z dehned below to data samples where / and / are not antiparallel. 

We hrst discuss the case of a Higgs boson heavy enough to decay into ti pairs. We 
dehne a sample A containing events where the top quark decays semileptonically and the 
top antiquark decays hadronically, 




t ^ + b 

i — W~ + b^q + q' + b. 


(HI.l) 


The sample A is dehned by the charge conjugated decay channels of the tt pairs, 




t —>• W~^ + b^q + q' + b 
t —^ Pfi b . 


(III.2) 


Each of these samples contains a fraction of 2/9 of all ti pairs. If we dismiss r leptons as 
spin analyzers for the top quark, the remaining fraction is 4/27. 

The above decay modes are especially suited to study correlations that result from top 
spin-momentum correlations. From the hadronic decays of the t (i) the momentum of the 
t {i) may be reconstructed on an event by event basis and the acompanying lepton of the i {t) 
decay serves as an excellent spin analyzer. The knowledge of the top momentum allows a 
reconstruction of the top rest system. In our approximation we have kj = —k^, which means 
that also the t and i rest system is known in the case of (|HI.1|) and ([HI.2|) , respectively. We 
may therefore dehne the CP-odd correlation 


■ p;+ y + (ki • p;- y. (in.3) 

where p^+ is the hight direction of £~^ in the top quark rest system, and p|_ is the hight 
direction of £~ in the top antiquark rest system. (In the more general case one should 
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replace p|+, p|_ —> p£+, pi- defined in the ip rest frame.) To measure £i, the corresponding 
expectation values are taken separately with respect to the two samples A and A. The 
correlation (|111.3|) traces the CP-odd spin asymmetry ( C>i). 

We further dehne the following triple correlation: 


^2 = (ki ■ (p;+ X Pe) - (ki ■ (Pf- X pj) )^ , 


(III.4) 


where p^ is the unit momentum of the b in the i rest system and p^ is measured in the t rest 
system. The triple correlation (|III.4|) probes the CP violating spin-spin correlation { 02 ). 

The calculation of the correlations ( [tll.3|) , ([tll.4|) involves a trace over the spin spaces 
of t and t of the form Tr{i?pi ® pi}. The decay density matrix pt for semileptonic t decays, 
which is needed for expectation values taken with respect to sample A, is given by 


pt{t = 


6x+(l— x+) 


1 - 3^2 + 2^3 


11 -h P^+ ■ (Tt 


dx. 


dVLi+ 

dvr 


(III.5) 


where k = rn^/m^ and x+ = 2E^+/mt G [k, 1] is the scaled energy of the lepton. Further 
we need the decay density matrix for hadronic t decays, in which the b analyzes the top 
antiquark spin. It reads: 


Pi{t -^W b ^ q + q' + b) = 


fl H-Ph ■ CTf 


dVLi 

dvr 


(III.6) 


The corresponding decay density matrices for sample A can be obtained by replacing x+ —> 
x_, p}+ -p|-, dVL(,+ dQe-, crt crt in (|III.5|) , and replacing pg ^ -p^, CTf cr*. 


dVti —> dklfc, (Tt (Tt in ([1II.6I) . Using these decay matrices, we get for the correlations £ 1 ^ 2 '- 


£1 = 
£2 = 



8 

9 


I-2k 
1 + 2k 


( 02 ) . 


(III.7) 


The factor (1 — 2 k)/(1 -|- 2k) ~ 0.41 is the spin analyzer quality of the b{b). The statistical 
sensitivities of these correlations on the unknown couplings at and dt are determined by the 
mean square fluctuations of the observables. We have 
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(III.8) 


((kf 


((ki ■ P^)^ ^ ’ 

((k,-(p|-xp*))^)_^ 


2 

9 ■ 


In order to decide whether or not the Higgs boson is a CP eigenstate, it is sufficient to 
observe a nonzero value of or £ 2 , which implies a nonzero (Oi) or ((P 2 ), respectively. 

fl 2^ — 

The number ' of (p ^ tt events that are needed to establish a nonzero correlation 
( ) with three standard deviation (s.d.) significance are given by: 



= 9- 

27 

3 

1 

4 

2 

(o,y ’ 


= 9- 

27 

9 

/I + 2k 

4 

16 

U-2k 


02 ^ 


(III.9) 


Here we took into account only electrons and muons in the samples A and A. In Fig. 10 we 
plot these numbers as a function of r* defined in ( |H.17| ) for = 400 GeV (again assuming 
without loss of generality atdt > 0). For a given number of events we can read off from the 
figure the interval for which leads to a nonzero correlation with 3 s.d. signihcance. For 
example, values 0.18 < rj < 0.52 would generate a nonzero £2 with 3 s.d. signihcance in a 
data sample of = 1500 events. If no effect is seen in this sample, this interval for is 
excluded. 

The above CP studies may be complemented by considering the correlation 


^3 = ( P^ ■ Ps )_4 + ( Pe- ■ Pb ) j , 


A 


which is related to the CP-even spin-spin correlation (O 3 } by 

1 -2k 


& = 


9 1 + 2k 


(O 3 } . 


From ( (p*+ ■ P5)2 = ( (p*_ ■ ) _ = 1/3 we get 


A 


“ 4 32 U-2Ky (Os)^ ■ 


(HI.IO) 


(HI.ll) 


(HI.12) 
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This number is also plotted (dotted line) as a function of r* in Fig. 10. A simultaneous 
measurement of £2 and £^ with Nti = 1500 would have a 3 s.d. sensitivity to values of 
between 0.18 and 1 . 


Higgs boson production processes where the tp rest system can be reconstructed include 


e+e“ —>• Zip and e+e" —^ e^e~p at a future linear collider and, if realizable, —> 

p. For the W^W~ fusion process —>■ and for p production at the LHC the 

above observables should be used with momenta obtained in the laboratory frame P|. This 
decreases the sensitivity of the observables as compared to above. 

Essentially the same analysis as above can be carried out for the case p —> . This 

mode is of interest for a quantum number analysis of Higgs bosons with mass < 2mw 
For a wide range of model parameters one typically expects in this case the branching ratio 
of the p ^ T~T~^ mode to be about 8 percent. For analyzing the r spin several decay modes 
can be used. We discuss here only the decays and p^z/.r(W)- The 

decay density matrices are of the form 

d£lg^ 




11 ± CbPI^ ■ (TrT 


dvr 


{B = 7r,p) , 


(HI.13) 


with Cjr = 1, Cp = 0.456 and is the direction of flight of B^ in the rest system. 
The mode r —>• Oi, respectively t ^ Sn can also be taken into account, see [^. We will 


assume in the following that the flight directions can be reconstructed for the above 
decays. This should be possible with some effort, as was done in the CP studies in 
production at LEP p3|-|35| (see also [^). Since the CP-odd spin asymmetry (Oi) is tiny 
for p —>■ we only discuss the CP-violating spin-spin correlation ( O 2 ) and the CP-even 

quantity ((P 3 ). Table 2 summarizes our results for the different decay modes of the 
pairs and the corresponding correlations of the final state momenta which are sensitive to 
{O 2 ) and {Os). 

To constrain the unknown parameter Vr, we combine all these decay channels. The 
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number of (p ^ t events to get a 3 s.d. significance for a nonzero ( O 2 ) is: 

1 1 


+ = 9--■ 


(III.14) 


8 (C>2)2 0.013 + 0.056 c2 +0.064 c4 ■ 

is obtained from (p.ll.l4|) by the replacements 9/8 —>• 27/16 and {O 2 ) {^ 3 )- 

Both numbers are plotted in Fig. 11 as a function of r,-. According to Fig. 11 about 2000 
ip —> T~T~^ events would be needed to establish r^- = 0.^1 [a^ + 0 ^-) > 0.36 also at the 3a 

level. 


IV. CONCLUSIONS 

The analysis of the decays p ^ ff made above shows that the spin-spin correlations 
O 2 and C> 3 , respectively the correlations £2 and £3 are useful tools for determining the CP 
properties of neutral Higgs bosons. If there is CP violation in the Higgs sector it could 
be established with these observables in a direct way. For a light Higgs boson with mass 

< 2 mw the decay p ^ T~r^ offers a good possibility to determine whether or not p is 
a CP eigenstate. We found that with 2000 (p —> r“r+ events the combined use of £2 and £3 
yields a sensitivity at the 3 ct level to a range of scalar and pseudoscalar Yukawa couplings 
corresponding to the ratio = 0.^1 {pr + o^t) > 0.36. The sensititvity can be improved by 
using also other r decay modes than those considered above. 

For the case of p ^ ti we showed that the order ag QCD corrections to the single 
spin asymmetry and to the spin-spin correlations are small. Again, a sizeable pseudoscalar 
component of p can be traced with the combined use of the correlations £2 and £ 3 . We 
found that r* = at/{at + at) > 0.18 could be established as a 3a effect with about 1500 
reconstructed p —>■ ti events. 

If Higgs boson(s) will be found these quantum number analyses should eventually be 
feasible with appropriately tuned new high luminosity colliders. 
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Table 1: Transformation properties of the 
strnctnre fnnctions nnder discrete 
symmetry transformations. 



Table 2: Different decay modes of the r pairs and the corresponding correlations of 
the hnal state momenta which are sensitive to {O 2 ) and {O 3 ). 
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Figure la: Decay rate ip ^ tt in GeV as a function of the Higgs mass 
for mt = 175 GeV, at = 1 and at = 0. The dashed line 
represents the Born result and the full line is the result to 
order a<,. 




Figure lb: Decay rate (p ^ ti ui GeV as a function of the Higgs mass 
for mt = 175 GeV, at = (2/3)^/^ and at = (1/3)^/^. The 
dashed line represents the Born result and the full line is the 
result to order Ug- 




Figure Ic: Decay rate ip ^ tt in GeV as a function of the Higgs mass 
for mt = 175 GeV, at = 0 and at = 1- The dashed line 
represents the Born result and the full line is the result to 
order a„. 





GeV 


Figure 2: Expectation value ( Oi ) as a function of the Higgs mass for 
mt = 175 GeV, at = (2/3)^/^ and dt = (1/3)^/^. The dashed 
line represents the Born result and the full line is the result 
to order 
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Figure 3: Expectation value ( Oi ) as a function of rt = at/{at + at) 
for fixed Higgs mass m^p = 400 GeV and mt = 175 GeV. 
The dashed line represents the Born result and the full line 
is the result to order a^. 




GeV 

Figure 4: Expectation value ( O 2 ) as a function of the Higgs mass for 
mt = 175 GeV, at = (2/3)^/^ and dt = (1/3)^/^. The dashed 
line represents the Born result and the full line is the result 
to order 
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Figure 5: Expectation value (O 2 ) as a function of n = at/{at + at) 
for fixed Higgs mass m^p = 400 GeV and mt = 175 GeV. 
The dashed line represents the Born result and the full line 
is the result to order 




0.25 



Figure 6: Expectation value (O 3 ) as a function of rt = at/{at + at) 
for fixed Higgs mass m^p = 400 GeV and mt = 175 GeV. 
The dashed line represents the Born result and the full line 
is the result to order as- 
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Figure 7: Expectation value (O 4 ) as a function of rt = at/{at + at) 
for fixed Higgs mass m^p = 400 GeV and mt = 175 GeV. 
The dashed line represents the Born result and the full line 
is the result to order as- 





Figure 8: Expectation value ( O 2 ) for (p as a function of 

Vr = dr/{ttr + Or) for fixed Higgs mass = 100 GeV. The 
dashed line represents the Born result and the full line is the 
result to order a. 




Figure 9: Expectation value ( O 3 ) for (p as a function of 

Vr = dr I {ttr + Or) for fixed Higgs mass = 100 GeV. The 
dashed line represents the Born result and the full line is the 
result to order a. 





Figure 10: Number of events ip ^ ti lo establish a nonzero correlation 
(C>i, 2 , 3 ) (with 3 s.d. signihcance) as a function of rt = 
at/{at + at) for hxed Higgs mass rtpp = 400 GeV and mt = 

175 GeV. The dashed line represents the result for the 

( 2 ) 

full line is the result for NT and the dotted line is the result 
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Figure 11: Number of events ip to establish a nonzero cor¬ 

relation (£> 2 , 3 ) (with 3 s.d. signihcance) as a function of 

r-r = ar/{aT + hr) for hxed Higgs mass m^p = 100 GeV. The 

( 2 ) 

full line is the result for N _ , and the dotted line is the 

r 

result for N^'^} +. 

T 
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